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Dust associated with various stellar sources in galaxies at all cosmic 
epochs remains a controversial topic, particularly whether supernovae 
play an important role in dust production. We report evidence of 

dust formation in the cold, dense shell behind the ejecta-circumstellar 
medium (CSM) interaction in the Type la-CSM supernova (SN) 2018evt three 
years after the explosion, characterized by a rise in mid-infrared emission 
accompanied by an accelerated decline in the optical radiation of the SN. 
Such a dust-formation picture is also corroborated by the concurrent 
evolution of the profiles of the Ha emission line. Our model suggests 
enhanced CSM dust concentration at increasing distances from the SN as 
compared to what can be expected from the density profile of the mass loss 
from a steady stellar wind. By the time of the last mid-infrared observations 
at day +1,041, a total amount of 1.2 + 0.2 x 10° M, of new dust has been 
formed by SN 2018evt, making SN 2018evt one of the most prolific 

dust factories among supernovae with evidence of dust formation. 

The unprecedented witness of the intense production procedure 

of dust may shed light on the perceptions of dust formation in 

cosmic history. 


The content and species of dust grains that are associated with stellar consideringambientgrainsinanyoutflowofstellarwindofthesupernova 
sources in galaxies at all cosmic epochs remain a controversial topic, (SN) progenitor may become immediately sublimated and destroyed by 
particularly whether supernovae play an important role in dust pro- the energetic radiation pulse produced by the SN explosion?*. To date, 
duction. Moreover, they may even carve dust-hostile environments!*?, freshly formed dust has been observed in a handful of core-collapse 
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Fig. 1 | Evidence of the presence of dust in SN 2018evt. a, MIR and B-band light 
curves (black dots) of SN 2018evt. All phases are given relative to the estimated 
B-band maximum at MJD 58352. The Spitzer and NEOWISE observations are 
shown by purple diamonds and red squares as labelled. The purple and red 
curves fit the MIR band 1 and band 2 photometry before and after day +310, 
separately. The black line fits the linearly fading B-band photometry before day 
~400, with a decline rate of 0.624 + 0.006 mag 100 day”. b, Red-to-blue EW ratios 
of Ha (green circles), PaB (red stars) and Bry (blue triangles) lines. c, Evolution of 


the flux-weighted centroid velocities AV of Ha, Paß and Bry lines labelled with the 
same symbols as b. The AV of Ha measured before and after day +310 are fitted 
with separate linear functions as displayed by the two green line segments. 

d, Evolution of the EW of the Ca 11 IR triplet (Extended Data Table 2) and Ha lines. 
For the purpose of the presentation, the EW of Ha has been multiplied by a factor 
of 2. The error bars shown represent lo uncertainties of magnitudes, EW ratio, 
centroid velocity and EW. 


(CC) supernovae, both in the ejecta in situ> * and its interaction zone 
with the circumstellar medium (CSM) (for example, refs. 5,9,10). 
No clear observational evidence thus far shows any major formation 
process of dust grains in the thermonuclear runaway of -1 M, carbon/ 
oxygen white dwarfs (WDs)™”. 

Type la supernovae are generally thought to result from ther- 
monuclear explosions of WDs in binary systems. A rare subclass of la 
supernovae is denoted SN Ia-CSM, which is thought to be an exploding 
WD surrounded by a substantial amount of CSM”. The spectra of such 
events near peak luminosity are characterized by narrow Balmer emis- 
sion lines superimposed together with relatively shallow Fe-group and 
intermediate-mass elements. SN 2002ic was the first reported case 
of ala-CSM SN that revealed large amounts of CSM seen as a strong 
hydrogen emission“ ™. A number of additional la-CSM supernovae 
have been discovered and studied in detail in recent years, which 
include supernovae 2005gj, PTF11kx, 2012ca, 2013dn and 2015cp and 
a recent sample of Zwicky Transient Facility supernovae (ref. 18 and 
references therein). 

SN 2018evt (ASASSN-18ro (ref. 19)) is a la-CSM SN found in the 
spiral galaxy MCG-01-35-011 at redshift z = 0.02523 (ref. 20). SN 


2018evt shares some common optical spectral features with typical 
Type la SN 1991T-like supernovae, as shown in Extended Data Fig. 1a. 
They are characterized by strong Fe 111144.04 and 45129 absorptions, 
visible Si 1114564, weak S 11 W and Si 1116355 and lacking absorption 
features of Ca 11 H and K and Ca I! infrared (IR) triplet before maxi- 
mum optical light”’”’. The early phase light curves of SN 2018evt are 
comparable to those of SN 1991T, as shown in Extended Data Fig. 1b. 
The power-law fit of the earliest light curve of SN 2018evt (<-10 days) 
suggests a rise time t, = 18.76 + 0.24 days, whichis consistent with that 
of SN 1991T/1999aa-like events”. The inset of Extended Data Fig. 1b 
shows the early phase B - Vcolour curve, whichis also in general agree- 
ment with that of SN 1991T after correcting the host reddening of 
E(B - V) < 0.32 mag, which has been estimated from the equivalent 
width (EW) of the NaI D lines”. The presence of the Ha line makes it 
ala-CSM SN similar to SN 2002ic. The entire spectral sequence of SN 
2018evt directly resembles other well-observed la-CSM supernovae 
events such as PTF11kx and SN 2002ic (Extended Data Fig. 2). The 
near-infrared (NIR) spectrum of SN 2018evt at ~324 days after the 
maximumiis similar to that of another la-CSM SN candidate SN 2012ca 
with data at a comparable epoch (Extended Data Fig. 3). 
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Fig. 2| SED fitting of SN 2018evt at the rest-frame wavelength. The optical- 
to-NIR BVgriJHK, data are fitted by a single BB before day +674. Emissions from 
the CSM dust calculated from our double-shell model (see ‘BB fit and dust 
sublimation’ and Fig. 4 for more details) are illustrated by cyan-dashed lines. 


Wavelength (um) 


Emissions from the newly formed dust are shown as yellow-dotted lines. Note 
that the thermal emission of the newly formed dust becomes progressively 
more dominant over time after day +445. The error bars shown represent 10 
uncertainties of the monochromatic luminosities. 


Results 
Observations. SN 2018evt shows the characteristic spectral features 
ofa Typela-CSM SN at early times together with Balmer emission lines 
(Extended Data Fig. 1) and the typical long-duration optical/IR light 
curves at late phases (Extended Data Fig. 4), indicating a continu- 
ous interaction between the expanding ejecta and a radially extend- 
ing CSM. We observed SN 2018evt with the Spitzer” InfraRed Array 
Camera (IRAC) at 3.6 um (CH1) and 4.5 um (CH2)” in 2019 (Extended 
Data Fig. 5a and Extended Data Table 1). Meanwhile, the area of the SN 
location was scanned by the Near-Earth Object Wide-field Infrared 
Survey Explorer (NEOWISE) reactivation mission” at 3.4 um (W1) and 
4.6 um (W2) from 2019 to 2021” (Extended Data Fig. 5b and Extended 
Data Table 1). The mid-infrared (MIR) fluxes of SN 2018evt exhibit an 
initial decline from +149 to +310 days relative to the estimated B-band 
maximum at modified Julian day (MJD) 58352; however, it is followed 
by anunprecedented rebrightening until the SN reached its peak lumi- 
nosity in both the W1 and W2 bands at around day +674 (Fig. 1a). This 
behaviour is not only distinct from the steadily fading light curves in 
optical bandpasses but also has not been seen in any previous la-CSM 
supernovae in similar MIR filters (Fig. 1 and Extended Data Fig. 4b,c); it 
is, however, likely that this is due to the lack of adequate time coverage 
of the observations of the latter. 

The optical spectral sequence of SN 2018evt spans days +125 
to +579 and also reveals conspicuous temporal evolution of the 


asymmetric characteristic Ha profile. We measure the EW separately 
for the red and the blue wings of the Ha (shown in Extended Data Fig. 2, 
Extended Data Table 2 and ref. 20), PaB and Bry profiles (Extended Data 
Fig. 3 and Extended Data Table 2). The ratios of red-to-blue wing flux 
increase steadily from day +125 to ~+310 but turn over and decrease 
afterward (Fig. 1b), in pace with the MIR flux evolution. Meanwhile, 
the flux-weighted centroid velocity AV of the Ha line (see ‘Analysis 
of the spectroscopic behaviours of SN 2018evt’ for details) evolves 
steadily from the blueshifted side (-400 km s”) to the redshifted side 
(+300 km s”) before day +310 and thereafter moves gradually back to 
the blue side (-200 km s”) (Fig. 1c). In addition, the evolution of the EW 
ofthe Ca I1 NIR triplet also exhibits a fall and rise, in concert with the evo- 
lution of the MIR flux and the Ha line profile (see, for example, Fig. 1d). 


Model. The slowly declining luminosity in the optical and NIR (Fig. 1 
and Extended Data Fig. 4) and the broad, long-lived Ha line (as shownin 
ref. 20 and Extended Data Figs. 2 and 6) dominating the late-time spec- 
tra of SN 2018evt both indicate that a substantial amount of late-time 
emission would arise from kinetic energy from the ejecta—CSM inter- 
action converted to radiation”. Such an additional energy source 
leads to a much slower luminosity decline (Fig. 1a) as powered by the 
5éCo > Fe decay: that is, -0.97 mag 100 day". In such a context, a cold, 
dense shell (CDS) develops during the ejecta—CSM interaction, with 
the CDS being located at a region between the shocked CSM and the 
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Fig. 3 | Time evolution of the different radii. The BB photospheric radius Rpg 
and the BB dust radius RN" are derived by fitting a BB spectrum to the optical-to- 
NIR luminosity and the MIR flux excesses, respectively. The latter is displayed in 
Fig. 4, and the associated temperature of the newly formed dust can be seen from 
the inset of Fig. 6. The horizontal grey-dashed line indicates the inner radius of 
the inner shell of the CSM (2.2 x 10° cm) in the double-shell model. The inset 
presents the temporal evolution of the FWHM width of the broad Ha line. Two 
blue line segments present linear fits of the data before and after day +310, 
respectively. The shock radius R, was derived by equation (1) in ‘BB fit and dust 
sublimation’ by assuming that the shock velocity was 10,000 kms ‘before the 
first observation at day -+120 and approximated by the FWHM width of the broad 
Ha afterward. The error bars shown represent lo uncertainties of FWHM. 


shocked ejecta**”’. This is the region where the SN ejecta and the CSM 
mix produce suitable conditions that allow the condensation of dust 
grains on short timescales’. 

Black-body (BB) fitting of the spectral energy distribution (SED) 
of SN 2018evt over the optical-to-NIR wavelength range suggests a 
broad temperature range of around 6,400-7,000 K during our observa- 
tions. As shown in Fig. 2, the MIR flux excess is obvious and becomes 


progressively more dominant over time as the optical emission 
decreases. The MIR flux excess can be attributed to the thermal emis- 


sion from dust at temperatures of 100-1,000 K (refs. 31,33). The pho- 
tospheric radius RODE estimated from the BB fitting is shown in Fig. 3. 


After day +141, ROP" decreases continuously with time (Fig. 3). This 
indicates a progressive deviation of the BB photosphere from the 
expanding CDS (for example, Fig. 7 of ref. 34), allowing the CDS to cool 
to a lower temperature. 

We explore various radial profiles of pre-existing CSM that may 
account for the time-variant excess of MIR flux due to the thermal emis- 
sion of dust. An initial decrease before day +310 could be attributed 
to a single-shell IR echo or a prominent process of dust sublimation 
as the forward shock runs through. The subsequent brightening after 
day +310 would suggest that newly formed dust accounts for the later 
MIR emission, in either the postshock regions of the CSM or the cooling 
ejecta. Assuming CSM dust density follows a power-law distribution 
Paus * r5, the plausible fit to the time-variant MIR flux excesses before 
day +310 (Fig. 4) by searching among a grid of parameters requires a 
power-law index s = 1.15. Other free parameters include the total optical 
depth and the inner and outer radii of the CSM shell (‘BB fit and dust 
sublimation’). The shallower radial density profile (s = 1.15) implies 
enhanced dust content at larger distances from the progenitor star. 

In the case of the steady dust mass loss s = 2, a plausible fit can 
also be achieved by introducing two shells of pre-existing CSM dust 
before day +310, namely the double-shell model. As the forward shock 


propagates outwards, grain sublimation takes place progressively only 
within the inner shell at a distance of 2.2 x 10*% cm, while the emitting 
dust grains in the outer shell, which is located at 6.0 x 10” cm from 
the SN, remain unaffected early on (Figs. 3 and 5). Because of the lack 
of early time spectral coverage, we adopt an initial shock velocity 
V,=10,000 km s” before day +120 based on the value estimated for 
SN 2002ic (refs. 16,35). As evidenced by the decreasing full-width at 
half-maximum (FWHM) width of the Ha profile (inset of Fig. 3), the 
forward shock expands into the inner shell of the CSM and deceler- 
ates. The progressive destruction of the inner shell dust grains leads 
to a continuously decreased emission in the MIR (Fig. 4). After -+310 
days, the forward shock supersedes the outer bound of the inner shell 
and enters a relatively low-density zone between the two CSM shells. 
The MIR emission becomes increasingly dominated by the relatively 
constant thermal emission from the outer shell. Our modelling sug- 
gests a massive outer shell of 5.2 x 10° M, of dust and an inner shell of 
3.2 x 10° M, of dust, corresponding to two episodes of elevated dust 
mass loss of 2.1 x 10° M, yr‘and1.8 x 107 M, yr”, respectively (‘BB fit 
and dust sublimation’). 

Akin to the single-shell model, our double-shell model also sug- 
gests enhanced dust concentration at larger distances from the SN as 
compared to what can be expected from the density profile of the mass 
loss froma steady stellar wind. The dust distribution inferred fromthe 
MIR flux excesses before day +310 can be modelled in terms of a double 
shell, which assumes a sudden change of the density profile of the dust, 
or a single shell with a flatter radial profile (Fig. 4). 

Both the single-shell and double-shell models are compatible with 
the MIR flux excesses at day $+310, but they cannot fit the MIR flux 
excesses at day >+310. After day +310, the rebrightening of SN 2018evt 
in the MIR demands notable contributions by additional emission 
sources, which can be well-attributed to the emergence of warm dust 
in regions behind the forward shock. As shown in Fig. 3, the BB radius 
Rip of the newly formed dust content fitted to the SED after day +310 
increases monotonically and remains within the shock radius R,. 

The inferred mass of the newly formed dust increases over time 
following a relation M,« t and reaches 1.2 + 0.2 x 10? M, by the 
time of our last observations at day +1,041 (Fig. 6 and Extended Data 
Table 1). The errors of the dust mass and temperatures are deduced 
using the Monte Carlo method via propagation of optical-to-NIR pho- 
tometric errors into BB fits and the MIR photometric errors into the 
flux-excess calculations. The dust sublimation timescale is extremely 
sensitive to the temperature close to the binding energy of the dust 
particles”. Dust survival close to the shock is possible if the dust dis- 
tribution is patchy or in an opaque disk, in which the self-shielding of 
the dust particles is important”. Our double-shell model assumes that 
a substantial amount of dust may survive the initial UV/optical emis- 
sion of the SN explosion out to the inferred inner CS dust shell radius 
of 2.2 x 10% cm, as shown in Fig. 3. 

Moreover, our model with dust formation is also consistent with 
the time evolution of the observed colours in the optical. The colours 
can be modelled by including the absorption and scattering effects of 
the newly formed dust (Extended Data Fig. 7). The increasing amount 
of dust after day +310 may contribute to the apparent blueward evolu- 
tion of the B-V, g-rand g-icolours. An increasing amount of scattered 
light is expected with more dust, which leads to excess flux inthe Band 
g bands, as shown in Extended Data Fig. 7. At even later epochs after 
day ~+500, the SN also exhibits accelerated fading in optical band- 
passes, which is compatible with a change from the optically thin to 
optically thick regimes of the newly formed dust. Such a transition 
is similar to the dust-formation process observed in the ejecta of SN 
1987A (ref. 38). 


Discussion 
Ha emission is powered by the interaction between the ejecta and the 
CSM""S, A thorough investigation of the time series of spectroscopy and 
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Fig. 4 | Modelling to the MIR flux excesses of SN 2018evt. The single-shell CSM 
dust model (solid green line) with an inner radius of 2.6 x 10” cm can fit well 

the declining flux excess in the MIR at day $+310 and infer a flatter power-law 
index of the dust density s = 1.15. In the case of the steady-wind mass loss s = 2.0 
for the double-shell model (solid black line), dust grains within the inner shell 

at 2.2 x 10% cm were continuously destroyed by the expanding forward shock 
between days -+200 and +310, causing a monotonically decreased flux excess 

in the MIR (dashed grey line). The presence of an outer CSM dust shell with an 
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inner radius of 6.0 x 10” cm is necessary to account for the time evolution of the 
flux excess before day +310 (dotted grey line). The prominent rise of the MIR flux 
excess of SN 2018evt after day +310, which cannot be explained by the thermal 
emission of any pre-existing dust content, demands a substantial amount of new 
dust to form promptly in the postshock regions (dotted red line). Panels (a) and 
(b) present the flux excesses of SN 2018evt at -3.5 um (Spitzer CH1 and NEOWISE 
WI) and ~4.6 um (Spitzer CH2 and NEOWISE W2), respectively. The error bars 
shown represent lo uncertainties of monochromatic luminosity excesses. 


spectropolarimetry within the first year of the SN explosion suggests 
that the SN ejecta expands into a dense torus of disk-like CSM”°. Such 
aconfiguration is in good agreement with the picture depicted by the 
spectroscopic and MIR flux evolutions that span days +125 to +1,041. 
The SN ejecta running into highly asymmetric disk-like CSM leads toa 
high-density torus inclined at an angle towards the observer. The early 
blueshift of the Ha emission line is explained if the redshifted side of 
the shocked CSM is blocked by the photosphere, as shown schemati- 
cally in Fig. 5. The redward shifts of the Ha emission line (that is, from 
days +125 to +310; Fig. 1) are caused by the receding photosphere as 
the photosphere shrinks, as proposed for the CSM configuration in 
PTFlliqb (ref. 39; their Figs. 10 and 12). After day +310, warm dust grains 
start to coagulate in the CDS and gradually block the receding side 
of the Ha line again, resulting in a blueward shift of the line profiles 
(Fig. Sc). After day +674, the W1-W2 colour of SN 2018evt becomes 
redder over time (Fig. 1a), indicating a decrease in the temperature 
and the MIR emission of the newly formed dust. 

The presence of a highly asymmetric ejecta-CSM interaction 
zone is also supported by detailed spectropolarimetry of SN 2018evt, 
which shows a wavelength-independent degree of polarization with 
non-evolving position angles that is characteristic of electron/dust 
scattering froma highly axisymmetric configuration®”°*°. Despite 
assuming spherical symmetry, both the single-shell and double-shell 
distributions of the CSM dust shell provide a satisfactory description 
of the SED evolution of SN 2018evt spanning days +149 to +310, in 
particular the time-variant excess in MIR. By incorporating the geo- 
metric information obtained from spectropolarimetry”’, the ejecta- 
CSM interaction process of SN 2018evt before day +310 is illustrated 
by the schematic sketches Fig. 5a,b. In our double-shell model, the 
dust in the CS wind at the radius ~2.2 x 10™ cm may be distributed in 
a disk or torus instead. The destruction and formation of the dust 


manifests qualitatively similar trends inthe temporal evolution of the 
MIR excess. 

The CSM masses derived from optical and optical-to-MIR luminosi- 
ties in shock interaction regions are -0.2-4.5 M, (see ‘The progenitor’s 
mass loss’ for details; ref. 20), corresponding to mass-loss rates of 
M =1x10-3to9x10°M, yr. Such CSM masses estimated from the 
kinetic-to-radiation energy process across the shock front appear to 
be -10° larger compared to the amount of dust within the inner CSM 
shell (3.2 x 10° M,), which contributes most of the MIR excess before 
day +310 (Fig. 4). Thus a very low dust-to-gas mass ratio within the inner 
shell at a relatively smaller distance (2.2 x 10° cm) can be inferred, 
which is likely caused by the prompt destruction of a substantial 
amount of grains inthe inner shell by energetic particles from the SN’. 
Onthe other hand, a gas-to-dust mass ratio on the order of 100 canbe 
inferred in the more massive (5.2 x 10° M,) and distant (6.0 x 10” cm) 
outer shell, whichis consistent with what is anticipated in the interstel- 
lar medium“. This probably means that the dust in the outer shell is 
much less affected by both the radiation field of the SN and the ener- 
getic particles from the shock interaction between the ejecta and the 
inner shell. Results similar to the outer shell can be derived by compar- 
ing the above CSM masses and the dust mass (6.0 x 10 * Mo) located at 
2.6 x 10” cm in the single-shell model. 

The mass loss of the progenitor before the explosion is in favour 
of either a thermonuclear explosion from a WD + asymptotic giant 
branch (AGB) star system’ or a core-degenerate system in which a 
WD merges with the core of a massive AGB star that triggers a thermo- 
nuclear explosion at the end of acommon envelope phase or shortly 
after”. The mass loss is also consistent with a WD + main sequence 
system for the common envelope wind model**”’. The progenitor 
systems are consistent with the measurements of the wind velocity 
V,,=91+58 kms “from the absorption minimum of the narrow P Cygni 
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Fig. 5 | Schematic sketches of SN 2018evt at different phases. The blue arrow 
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torus with pre-existing dust being destroyed by SN radiation and CSM shock. 
Geometric configuration before day +310 when the redshifted component of the 
Ha line (shown as solid blue ellipses) is blocked by the photosphere, producing 
blueshifted line profiles. b, As the SN ejecta expands and the photosphere 
recedes over time, more redshifted emission is revealed, resulting in a redward 
evolution of the line centroid as seen in Fig. 1c. c, When new dust forms at the 
postshock CDS (thick solid line), the redshifted side of Ha is blocked again, 
leading to blueward evolution of the line profile. 
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Fig. 6 | Temporal evolution of the mass of the newly formed dust. As shown by 
the black line, the mass of the newly formed dust of SN 2018evt can be well fitted 
by a power law: Ma ~ t* for 0.3 um graphite grains. The dust masses calculated 

for graphite and silicate particles of radius 0.05 um are also presented. The inset 


traces the temperature evolution of newly formed graphite dust particles of 
radius 0.3 um. The dust masses estimated for Type IIP supernovae 2004et and 
1987A and IIn supernovae 200Sip, 2006jd and 2010jl are also shown. The error 
bars shown represent lo uncertainties of masses and temperatures. 


profiles of the Ha line (Extended Data Fig. 6 and ref. 20). Compared 
with the density profile of the dust mass loss froma steady stellar wind 
s=2, a flatter radial profile s = 1.15 inferred in the single-shell model 
indicates enhanced dust concentration at increasing distances from 
the SN. The double-shell model also points to the same result that a 
higher dust mass-loss rate (M = 2.1 x 10-5 Mọ yr) of mass ejections 
is measured within the distant outer shell (6.0 x 10” cm) and a lower 
dust mass-loss rate (M = 1.8 x 1077 Mo yr~) is measured within the 
close inner shell (2.2 x 10% cm). Both the single-shell and double-shell 
models suggest an enhanced dust presence at larger distances from 
the progenitor star. This shallower radial density structure results from 
a variable mass loss, which is likely to happen in the entire AGB 
evolution’. 


The rebrightening in the MIR after day +310 can be modelled as a 
result of the formation of a substantial amount of warmer dust at late 
phases (Figs. 1, 2 and 4), distributed in a prolate shell vertical to the 
CSM disk (Fig. 5). It also provides a natural explanation of the 
red-to-blue emission-wing ratio of Ha due to uneven extinction by the 
newly formed dust (Fig. 1). This behaviour is also observed for the Ca 
1 NIR triplet, which can be similarly explained. The rapid weakening 
of the Ca 11 NIR triplet may also indicate the depletion of calcium by 
dust formation. The proposed process of dust formation is corrobo- 
rated by the time evolution of the EW of the Ca II NIR triplet emission 
lines, which exhibits an increase after day ~+310 (see, for example, 
Fig. 1d). Figure 6 shows the temporal evolution of the mass of the newly 
formed dust of SN la-CSM 2018evt in our model, compared with other 
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CC supernovae. The estimated dust mass is highly dependent on the 
species and size distribution of the dust grains. For graphite witha size 
of 0.3 um, the dust mass grows rapidly following a power law of index 
4 with the time after the explosion and reaches ~1.2 + 0.2 x 10°? M, at 
the last epoch of the MIR observations at day +1,041 (Fig. 6). For graph- 
ite or silicate of 0.05 um, the dust mass is about three or five times 
higher than the value derived for the 0.3 um graphite dust (Extended 
Data Table 1), respectively. The temperature of the newly formed dust 
is presented in the inset of Fig. 6. A monotonically decreasing tempera- 
ture from 1,000 to 500 K between day +434 and day +1,041 is likely to 
be mostly affected by the expansion cooling of the CDS region between 
the forward and reverse shocks. It may also be regulated by various 
heating mechanisms, including radiative heating from the SN shock, 
collisional heating with the ambient warm gas and energy exchange 
between the gas and dust*’. As the SN ejecta expands and the 
dust-forming region in the CDS cools, dust grains may continue to 
coagulate. Depending onthe duration of the timescale that the ejecta 
expands into the CSM, even orders of magnitude higher dust content 
may be produced during such a process. A notable fraction of the 
unburned carbon in the ejecta, if not all, can be locked in the newly 
formed dust. As progressively deeper layers of the ejecta move into 
the CDS, we may also expect a massive amount of iron and silicate dust 
to forminsuchan environment. The James Webb Space Telescope can 
probe the signatures of such dust in the coming years. The estimated 
mass (M,) and the BB emission radii RMR of the newly formed dust 
masses are consistent with those seen in CC supernovae at similar 
phases (Figs. 3 and 6 and ref. 30), suggesting a rapid and efficient 
mechanism for dust production in these supernovae. 

Finally, we remark that <1% of CC supernovae occur in elliptical 
galaxies in the local universe”, so dust production in thermonuclear 
explosion supernovae la is a major channel of dust enrichment in 
early-type galaxies. Type la supernovae may also contribute to the 
dust budget in spiral galaxies*®. SN 2002ic was the first la-CSM SN ever 
discovered and has a dwarf elliptical host“. The weak Ha of the host 
galaxy of SN 2018evt also implies overall less active star formation”. 
Even though la-CSM is a rare subclass of thermonuclear supernovae, 
the unprecedented witness of such intense production of dust grains 
may shed light on the perceptions of dust formation in cosmic history 
(‘The dust contributions of host galaxies by la-CSM SN events’). 


Methods 

Observations 

Early-phase observations of SN 2018evt. The early-phase observa- 
tions of SN 2018evt were conducted with the dual-channel optical/ 
NIR camera ANDICam on the Cerro Tololo Inter-American Observa- 
tory 1.3 m telescope. Two epochs of BVRI-band photometry were 
obtained on 13 August 2018 and 17 August 2018 before the SN was 
too close to the Sun. ANDICam has an optical field of view (FOV) of 
6.3’ x 6.3’ (0.37” pixel”) and a NIR FOV of 2.4’ x 2.4’ (0.27” pixel”). 
The extraction of the NIR-band photometry was not successful due 
to the lack of bright stars for astrometric calibration to combine 
the dithered images. BVR/ point-spread function (PSF) photometry 
was performed on the optical images using PSFEx”’ following the 
detailed prescriptions described by ref. 51. The PSF photometry was 
calibrated to the standard Pan-STARRS catalogue” of the bright- 
est field star at (RA, dec.) = (206.655661°, —9.680946°) (J2000) with 
B=16.495 + 0.034 mag, V= 15.812 + 0.012 mag, g = 16.031 + 0.002 mag, 
r=15.603 + 0.002 mag and i= 15.461 + 0.003 mag. The r- and i-band 
photometry of this field star has been converted to the standard 
Johnson R/system™ following the transforming equations provided 
by refs. 55-57. The early optical light curves of SN 2018evt are shown 
in Extended Data Fig. 1b. We also retrieve early time photometry of SN 
2018evt using the ATLAS**”’ forced photometry service in the cand o 
bands and the All-Sky Automated Survey for Supernovae (ASAS-SN°”) 
sky patrol interface. The background flux of the ASAS-SN data of SN 


2018evt has been estimated by the pre-explosion median flux recorded 
with the same aperture as used for the SN photometry, based ona total 
of 266 visits. Both ATLAS and ASAS-SN photometry are also shown in 
Extended Data Fig. 1b. 


Spitzer observations. SN 2018evt was observed (primary investi- 
gator (PI): Sijie Chen) with the Spitzer” IRAC at 3.6 um (CH1) and 
4.5 um (CH2)* at days +271, +286, +434 and +445. We utilized the 
level-2 post-basic calibrated data images from the Spitzer Heritage 
Archive, which were reduced by the Spitzer pipeline and resampled 
onto 0.6” pixels. Source detection and aperture photometry were 
performed on the images in Extended Data Fig. 5a without host sub- 
traction using SExtractor®. We remark that flux difference is less than 
or similar to 10% for the Spitzer/IRAC photometry with and without 
template subtraction in ref. 31, which is well within the photometry 
uncertainty. We applied aperture corrections following the IRAC Data 
Handbook. The level-2 post-basic calibrated data images have been 
calibrated in an absolute surface-brightness unit of MJy sr‘, which 
can be transformed into units of uJy pixel by a conversion factor 
of 8.4616 for the angular resolution of our IRAC images, 0.6” pixels. 
The flux was converted to absolute (AB) magnitude according to the 
definition m,, = -2.5 log10(f) + 8.9, where fis in units of Jy (ref. 63). 
The AB magnitudes of SN 2018 evt in the CH1 and CH2 bands are listed 
in Extended Data Table 1. 


NEOWISE observations. The SN 2018evt field was also observed by the 
NEOWISE reactivation mission in the W1 (3.4 um) and W2 (4.6 um) bands 
since late 2013 as an extension of the WISE ALL-Sky Survey”. Using 
the online version of the NEOWISE Image Co-addition with Optional 
Resolution Enhancement (ICORE)°*°°, we retrieve the co-added 
NEOWISE images that centred at SN 2018evt, with a FOV of 0.6° x 0.6° 
and resampled to a pixel size of 1.0”. Given that SN 2018evt explodedin 
August 2018, we take the co-added image from January 2017 to January 
2018 as the reference image for background subtraction and generate 
the difference images for every single-visit co-added image using the 
Saccadic Fast Fourier Transformation (SFFT)”. Extended Data Fig. 5b 
shows the NEOWISE reference and difference images at the position 
centred on SN 2018evt. The time series of the differenced images 
clearly shows the notable variations in the brightness of SN 2018evt. 
The signal was notable in January 2019. After a noticeable dimming in 
the next six months, a dramatic rebrightening followed in 2020. 

Aperture photometry was performed on the differenced images 
using SExtractor” and calibrated to the profile-fit magnitudes in the 
Vega system released in the ALLWISE Source Catalogue. The photo- 
metric errors were measured on the corresponding variance images 
and corrected by a factor of 2.75, which gives the ratio of the input to 
output pixel scale (Section 13 of ref. 66). Such estimated photometric 
error o for each visit is used if it is larger than the photon noise from 
direct photometry on the differenced images. The Vega magnitude 
of SN 2018evt was transformed into AB magnitude according to the 
magnitude offsets between the two systems®. The AB magnitudes of 
SN 2018evt in the W1 and W2 bands are listed in Extended Data Table 1. 
The MIR-band light curves of SN 2018evt are shown in Extended Data 
Fig. 4b,c, together with other la-CSM supernovae including 2002ic, 
2005gj (ref. 35), PTF11kx (ref. 69), 2012ca (ref. 70), 2013dn (refs. 31,71), 
and 2020eyj (ref. 72). 


Optical photometry at Las Cumbres Observatory. Extensive BVgri 
photometry spanning days +124 to +664 was obtained with the Sinistro 
cameras on the Las Cumbres Observatory 1 m telescope, a global net- 
work for SN observations. Images were bias subtracted and flat-field 
corrected using the BANZAI automatic pipeline. The background 
template was then subtracted from the preprocessed images, adopting 
the SFFT algorithm”. Finally, PSF photometry has been performed on 
differenced images using ALLFRAME”. We remark that the light curves 
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of SN 2018evt before day +365 were achieved without subtracting any 
background template and reported in ref. 20. Comparisons between 
their direct photometry and our template-subtracted photometry 
obtained at similar phases suggest good agreement. In particular, 
the systematic magnitude differences in BVgri yield -0.03 + 0.05, 
-0.03 + 0.05, 0.02 + 0.05, -0.06 + 0.04 and -0.14 + 0.06, respectively. 

Templates for our gri-band exposures were directly obtained using 
Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) 
cutout images for the gri bands. B- and V-band templates were con- 
structed using Pan-STARRS grimages with the formulaB=g + w x (g-r), 
V=g-wx (g-r). The parameter w was achieved by minimizing the 
global residual flux computed based onall field stars, which produces 
the cleanest subtraction. Thus the best coefficients obtained are w = 0.3 
for Band w= 0.5 for V, respectively. 

Zero-point calibration was conducted using local field stars by 
calculating the 30 clipped median of the differences between instru- 
mental magnitudes and the standard Pan-STARRS Catalog” for gri 
and the American Association of Variable Star Observers Photometric 
All-Sky Survey data release 9 catalogue” for BV (using only stars with 
magnitude within 10-18 mag and photometric errors o < 0.1 mag). 
The BV gri-band light curves of SN 2018evt are shown in Extended 
Data Fig. 4a. 


Optical photometry with XLST and LJT. Optical photometric observa- 
tions of SN 2018evt were also conducted with the 60/90-cm XingLong 
Schmidt telescope (XLST) of the National Astronomical Observatories 
of China (NAOC) under a long-term Tsinghua University-NAOC Tran- 
sient Survey” and the Yunnan Faint Object Spectrograph and Camera 
(YFOSC)” mounted on the 2.4-m LiJiang telescope (LJT) at the Yunnan 
Astronomical Observatories. SN 2018evt was observed in the imaging 
mode of YFOSC. Images obtained by the XLST and LJT were processed 
using an automatic custom pipeline based on the Image Reduction and 
Analysis Facility (IRAF). The pipeline reduction follows standard proce- 
dures including bias and flat-field corrections, astrometric registration, 
template subtraction and PSF photometry. The BVgri photometry is 
also shown in Extended Data Fig. 4a. 


Optical spectroscopy. We also obtained 17 optical spectra of SN 
2018evt. A log of the spectroscopic observations is presented in 
Extended Data Table 2. Six spectra were taken with the 3.6 m ESO Faint 
Object Spectrograph and Camera v.2 (ref. 77) mounted on New Tech- 
nology Telescope at La Silla Observatory during the extended-Public 
ESO Spectroscopic Survey for Transient Objects (ePESSTO)’*. The 
observations were carried out under ESO programmes 199.D-0143 
(PI: Smartt) and 1103.D-0328, 106.216C (PI: Inserra). Four spectra were 
taken with the YFOSC/LJT in long-slit spectroscopic mode, and three 
were taken with the Beijing Faint Object Spectrograph and Camera”? 
mounted onthe 2.16 m Xinglong telescope. One spectrum was obtained 
with the Wide Field Spectrograph (WiFeS) mounted on the Australian 
National University 2.3 m telescope at Siding Spring Observatory®°. 
Three additional spectra obtained at days +490, +516 and +531 were 
acquired with the Folded Low Order whYte-pupil Double-dispersed 
Spectrograph (FLOYDS”) mounted on the 2.0 mtelescope at Las Cum- 
bres Observatory (Extended Data Table 2). The twin robotic FLOYDS 
spectrographs are mounted on the Faulkes Telescope South at Siding 
Spring Observatory and on the Faulkes Telescope North at Haleakala. 
Apart from the three late-time spectra mentioned above, another 12 
spectra (Extended Data Table 2) were obtained with the same telescope 
from days +125 to +365 (ref. 20) and were also included in this paper to 
measure the Ha line profile (for example, red-to-blue emission-wing 
ratio and flux-weighted centroid velocity AV). The photometry of SN 
2018evt obtained with the global network of 1 m telescopes and the 
12-epoch FLOYDS spectroscopy before day +365 have been published 
in ref. 20, which focuses on the early ejecta-CSM interaction and the 
spectropolarimetric properties of SN 2018evt. 


All optical spectra were reduced using standard IRAF routines. Flux 
calibration of the spectra was carried out using spectrophotometric 
standard stars observed at similar airmass on the same night. The 
spectra were further corrected for atmospheric extinction using the 
extinction curves of local observatories. 


NIR spectroscopy. This Article includes eight NIR spectra (Extended 
Data Fig. 3 and Extended Data Table 2). Four NIR spectra were obtained 
with the medium-resolution 0.8-5.5 um spectrograph and imager 
on the 3.0 m NASA Infrared Telescope Facility on Mauna Kea, named 
SpeX®”. Two NIR spectra were acquired with the Folded port Infra- 
Red Echellette (FIRE) spectrograph® on the 6.5 m Magellan Baade 
telescope. Another two spectra were obtained with the Gemini NIR 
spectrograph (GNIRS)* on the 8.2 m Gemini North telescope. The 
SpeX, FIRE and GNIRS spectra were reduced with Interactive Data 
Language codes, Spextool®, firehose® and the XDGNIRS pipeline®*”’, 
respectively. 


Analysis of the spectroscopic behaviours of SN 2018evt. All spectra 
were corrected for the redshift z= 0.02523 of the host galaxy” and 
extinction from the Milky Way E(B - V) = 0.05 mag (ref. 88). Three 
spectra lines were normalized with a pseudocontinuum by linear fit- 
ting to the spectra ranges [6250, 6350] A and [6700, 6800] A for Ha, 
[12200, 12500] A and [13200, 13500] A for Pa and [21100, 21400] A 
and [21900, 22200] A for Bry. Thus Ha, PaB and Bry are located at 
[6350, 6700] A, [12500, 13200] A and [21400, 21900] A, respectively. 
All spectra were scaled to match the photometry in the optical band- 
passes at corresponding phases and further used to measure the Ha 
luminosity (Extended Data Table 2) and EW (Fig. 1d). For each flux 
spectrum, following an approach similar to the analysis of ref. 20), we 
fita double-component Gaussian function to the Ha profile to decom- 
pose it into a broad and an intermediate component. We found that the 
centre of the intermediate Gaussian component, which has a typical 
FWHM width of -2,000 km s”, shows only moderate shift over time until 
the last epoch of spectroscopy at day +579. Such behaviour is in overall 
good agreement with the analysis based on the spectra obtained before 
day +365 by ref. 20. The determination of the centre of the intermediate 
Gaussian componentalso allows us to compare the blue and red wings 
of several major emission lines. In particular, we present the red-to-blue 
EW ratios measured between the red and the blue wings for Ha, PaB 
and Bry features in Fig. 1b. Note that the determination of the centre 
of the Paß and Bry lines was carried out based ona single-component 
Gaussian fit due to the relatively low signal-to-noise ratio (S/N). The 
EW ratios for the Ha and Paf lines were computed over a velocity range 
of -8,000 to +8,000 kms“. A narrower velocity range of -3,500 to 
+3,500 kms ‘was used for the measurement of the Bry profile. 

After correcting for the redshift of the host galaxy, we define the 


Apeak—Ao 


flux-weighted centroid velocity as AV = a XC, where c gives the 


speed of light and A, represents the rest wavelength of the line centre: 
that is, Ag(Ha) = 6563 A. The flux-weighted peak wavelength, A,eax, is 
calculated as Apeak = r where / denotes the wavelength of any spec- 


tral element over the emission profile. Such a quantity weights each 
spectral element by its flux f, thus providing a more robust trace of the 
bulk velocity of the line-emitting zone. Figure 1c shows the AV derived 
for the Ha, PaB and Bry lines for SN 2018evt. In Fig. 1, the uncertainties 
of the EW ratio and centroid velocity were calculated through the 
Monte Carlo method, assuming that all spectra have 10% flux 
uncertainty. 

The day +307 WiFeS® spectrum obtained with a higher spectral 
resolution (R = 3000) presents a well-resolved narrow Ha P Cygni pro- 
file (Extended Data Fig. 6). A two-component Gaussian fitting pro- 
cess suggests that the FWHM widths of the broad and intermediate 
components are 5,877 + 32 kms ‘and 1,643 + 12 km s”, respectively. 
After subtracting the broad and the intermediate components from 
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the day +307 WiFeS spectrum, we fit the residual spectrum with two 
separate Gaussian functions to better separate the narrow absorption 
and emission components of the P Cygni profile. We inferred a redshift 
z= 0.02561 + 0.00019 by assuming the narrow-emission component 
peaks at the rest wavelength of Ha (inset of Extended Data Fig. 6). The 
wind velocity, which is measured from the blueshifted absorption 
minimum, gives V, = 91+58 km s”. Our measurements are consistent 
with those reported by ref. 20 within the uncertainties: for example, 
z= 0.02523 + 0.00015 and V, = 63 +17 km s”. The redshift values derived 
in both studies are also consistent with those reported to the NASA/ 
IPAC Extragalactic Database®’”’. Therefore, we used z = 0.02523 and V, = 
63km s” (ref. 20) throughout the Article due to the smaller uncertainty. 


BB fit and dust sublimation 

The effective BB temperatures and radii ROP' were estimated by fitting 
a BB curve to a time series of SEDs constructed from the optical (BVgri) 
and/or NIR (/HK,) light curves of SN 2018evt. Optical photometry was 
obtained by the global network of the 1 m telescope at Las Cumbres 
Observatory, and NIR photometry was taken from ref. 20. The latter 
spans days +141 to +314 and was obtained with the Gamma-Ray burst 
Optical/Near-Infrared Detector”! mounted on the 2.2 m Max Planck 
Gesellschaft/European Southern Observatory telescope, operated at 
the La Silla Observatory in Chile. All SEDs were constructed after cor- 
recting for the E(B - V) = 0.05 mag Galactic extinction®®. We adopt a 
distance of 103.3 Mpc for SN 2018evt following the rationale provided 
inref. 20. Owing tothe lack of early time data before day +141, we adopt 
the optical and NIR light curves of the well-sampled Type Ia-CSM SN 
2005gj (ref. 92) to generate the SED of SN 2018evt during the missing 
phases and thereafter to calculate the MIR emission of CSM dust 
through the absorption and re-emission processes. Such an approxima- 
tionis validated by the high similarity in prepeak and around day +140 
spectra, light curve shapes and absolute brightness in optical and NIR 
bandpasses between supernovae 2018evt and 2005gj. In detail, super- 
novae 2018evt and 2005gj are both Type Ia-CSM objects (Extended 
Data Fig. land Figure 7 in ref. 93). They share a similar BB temperature 
and radius at day ~+140 based on their multiband photometry, as shown 
in Fig. 3 and Table 8 of ref. 92. Also, the peak fluxes of SN 2005gj are 
comparable with those of SN1991T (Fig. 7 of ref. 92), whose early phase 
spectrum and light curves match well with SN 2018evt (ref. 94), as 
shown in Extended Data Fig. 1 (see refs. 95-98 for SN1991T). 

For supernovae whose late-time emission is mostly dominated by 
strong ejecta—CSM interaction, their effective ROP' isexpected to coin- 
cide with the radius ofa thin CDS**”” located between the shocked CSM 
and the shocked ejecta. The RGP" of SN 2006gy reaches its maximum 
value at day -+115 (Fig. 7 of ref. 34), while the expansion of its CDS con- 
tinues as indicated by a rather constant FWHM width of the Ha, which 
traces the expansion velocity of the CDS. Similar behaviour is also seen 
inthe Type IIn SN 2006tf (Table 3 and Fig. 15 of ref. 100). Therefore, we 
suggest that the BB-fitted Re does nottrace the emitting radius of the 
CDS Reps. The latter can be represented by introducing a dilution factor 
Ç, which cannot exceed unity and decreases over time*. The true emit- 
ting radius of the CDS is given by Rcps= R99" /vZ (refs. 34,100). 

Figure 3 shows the Roe of SN 2018evt through BB fit to the 
optical-to-NIR photometry. At day +141, we measure gS =3x105cm. 
After day +141, Ree is decreasing in Fig. 3, indicating that Reps has 
departed from the corresponding BB radius at day +141 (<1). We 
adopted approximately Re = 3x 105 cm as the lower limit of the 
expanding CDS radius at day +141. 

For SN ejecta whose radial density profile follows an inverse 
power-law distribution, Pejecta* r”, the shock radius is given by 
equation (1) of ref. 101: 


R, =947 x 1080 —3 ( Vs )( t 


(n-3)/(n-2) 
(n — 3) \3,000 kms-1 en) 


cm (1) 


In the literature, the shock velocity V, is often approximated by 
the velocity corresponding to the FWHM width of the Ha emission 
line’©? 1°, Assuming a typical shock velocity of 5,000-10,000 kms? 
(refs. 16,105,106), which is consistent with the Ha velocity width meas- 
ured during our spectroscopic observing campaign on SN 2018evt 
between days +125 and +546 (inset of Fig. 3). By adopting ann = 8.5 
ejecta density profile estimated for SN 2002ic (ref. 17), we estimate a 
forward shock radius R, = 2.8 x 10° cm at day +310. 

This shock radius R, is less than the dust evaporation radii of 
~4-9 x 10"° cm for silicate or graphite dust'”, assuming the SN luminosity 
to be Lpa = 10® ergs“, which falls between the maximum for SN 2018evt 
(10° ergs”) and the peak luminosity of SN 2005gj (10**” ergs ‘in Table 8 
of ref. 92). This suggests that the pre-existing CSM dust in the single-shell 
modelat 2.6 x 10” cmand inthe outer shell of our double-shell model at 
adistance 6.0 x 10” cmare unlikely to be sublimated by the SN radiation 
as the dust temperature at 2.6 x 10” cm can only be heated to a tempera- 
ture of about 970 K, which is lower than the evaporation temperature 
of 1,500 K for silicate and 1,900 K for graphite!®®™!0, The shock radius 
R, is comparable to the inner radius of the inner CSM dust shell in our 
double-shell model (2.2 x 10° cm), indicating the dust grains within 
the inner CSM shell are likely to be destroyed by the forward shock, if 
they survived the initial pulse of the electromagnetic radiation of the 
SN explosion due toa patchy dust distribution or inan opaque disk”. 


MIR flux excess. The MIR flux excess compared to the best-fit BB SED 
for different epochs of observations is shown in Fig. 2. Two epochs of 
MIR observations were acquired with Spitzer CH1 (3.6 um) and CH2 
(4.5 um) at days +271 and +445. The MIR observations at days +286 
and +434 were not presented as they are nearly identical to the results 
for the day +271 and +445 observations, respectively (Extended Data 
Fig. 4 and Extended Data Table 1). Six more epochs of observations were 
acquired with NEOWISE W1 (3.4 um) and W2 (4.6 um) at days +149, +310, 
+517, +674, +881 and +1,041. 

MIR emission excess typically suggests the presence of warm 
dust. The MIR filters used to observe SN 2018evt provide rather com- 
plete wavelength coverage of observations spanning the peak of the 
thermal SED from dust with temperature spanning 100 < 7, < 1,000 K 
(refs. 33,108). The MIR emission excess has been explained by the forma- 
tion of new dust grains in a handful of CC supernovae, both inthe ejecta 
in situ’ ® and in the interactions between the ejecta and the CSM (for 
example, refs. 5,9,10). Alternatively, the MIR emission excess could origi- 
nate from the thermal IR emitted by dust particles that were presentin a 
CSM before the SN event. Such primordial dust grains may have formed 
in the expanding matter blown from red giant stars or AGB stars’. In 
addition to the thermal radiation of pre-existing CSM grains, our models 
also include emissions from any newly formed dust to account for the 
extreme MIR rebrightening of SN 2018evt after day +310. 


Modelling the emissions of SN 2018evt: a model with one or two 
primordial CSM shells and new dust formed in the CDS region. Dust 
particles inthe CSM absorb some of the UV/optical photons radiated 
during the explosion of the SN and its ejecta interaction with the CSM 
and re-emit the flux in the IR bands, producing an IR echo". Such an 
IR echo can be used to constrain the CSM dust properties around the 
SN, suchas their distribution, mass and composition. IR echo models 
for spherically symmetric CSM shells have been developed to account 
for thermal emission from pre-existing CSM dust, which provides a 
plausible explanation for the late-time excess in the observed IR light 
curves of Type la (refs. 115-117) and Type Il (refs. 110,114) supernovae. 
Thetime evolution of the IR echo is related to the ultraviolet and optical 
light curves of the supernovae. 

At any given time, a distant observer will see the IR echo located 
within an ellipsoid, with the SN and the observer lying at its two foci. 
Such an ellipsoid traces an iso-travel-time surface of the light emitted 
by the SN, which expands over time. The position of any point within 
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the ellipsoid can be expressed as (r, 0), where r denotes the distance 
from the point to the SN and 0 represents the scattering angle. For dust 
particles of radius a located at (r, 0) and an SN located at a distance D 
from the observer, the total flux emitted by the IR echo at time t gives 


2 Rout 
FOS 1 na(mB,|T4(r, 8, D1Q, Br 2) 
Rin 


where Rand R,,, are the radii of the inner and outer dust shell, respec- 
tively, n4 is the number density of the dust particles, and the Plank 
function B, at frequency v is determined by the dust temperature 
T,(r, 0, t), which can be estimated from the SN luminosity. Q, denotes 
the absorption and emission efficiency of the dust grains. The light 
curve of the emitting IR echo from dust distributed in a shell shows 
a plateau lasting for a period of 2R,,/c (where cis light speed) and fol- 
lowed by a decline for a period of time that is related to the radial extent 
of the shell. Such behaviour is similar to that reported for SN 2005ip 
(ref. 9) and several other supernovae that show strong ejecta-CSM 
interactions (for example, refs. 107,110,114,115). 

Assuming the CSM dust density around SN 2018evt follows an 
inverse power-law distribution paust% 7°, the MIR flux emitted by the 
CSM dust within a single shell can be derived from equation (2). The 
single-shell model is assumed to be spherically symmetric and described 
by four parameters: the inner and outer radii of the shell (Rin, Rou), 
the optical depth in the B band (Tg) and the power-law index (s). We 
initially set Rn in the range of 50 Id < Rn < 150 Id (where ld is a light day) 
to make sure the thermal radiation of the CSM dust declines between 
days +100 and +300. We run the single-shell model in tens of thousands 
of grid points based on the four parameters (Tz, Rin, Rwiatn, S) and obtain 
a group parameter (0.07, 100 Id, 80 Id, 1.15) to well fit the data (see the 
green lines in Fig. 4), where Ryiatn = Rou ~ Rin. A flatter radial profile of 
the CSM dust density was inferred from the single-shell model due to the 
smaller s compared with the value for the steady-wind mass loss of the 
progenitor system (s = 2). This suggests an increased dust concentration 
at increasing distances from the SN compared to what can be expected 
from the density profile of the mass loss froma steady stellar wind. The 
mass of the CSM dust within the single shell is derived to be 6.0 x 10° M,. 

In the case of the steady-wind mass loss s = 2, a plausible fit can 
also be achieved by introducing two shells of pre-existing CSM dust 
before day +310. The inner shell predicted by our double-shell model 
was caught by the forward shock at day ~+200 (see, for example, 
Fig. 3). Dust grains within the inner shell are thus gradually destroyed 
as the shock runs. The expansion velocity of the shock was assumed to 
be V,=10,000 km s” due to the lack of observations before day +120, 
followed by acontinuous deceleration as traced by the FWHM width of 
the broad Haas shown by the inset of Fig. 3. Before day +1,041, the outer 
shell of the CSM, which emerges at 6.0 x 10” cm, remains unaffected 
by the forward shock. Our double-shell CSM model also provides a 
satisfactory fit to the monotonically decreasing MIR flux curves before 
day +310. The SED fits to the BVgriJHK, and MIR-band photometry are 
illustrated in Fig. 2. The IR echo light curves of the two shells are shown 
in Fig. 4. The best-fit parameters are (Tg, Rin Rwiath) = (0.07, 8.5 Id, 2.2 Id) 
for the inner shell and (0.17, 230 Id, 30 Id) for the outer shell. The total 
optical depth of the pre-existing dust shells in the B band is Tg = 0.24, 
corresponding toa V-band extinction of Ay = 0.26 mag. We also remark 
that such an integrated extinction is consistent with the value estimated 
from the NaI D lines’. By assuming the CSM shells were built up by 
multiple epochs of pre-explosion eruptions, the derived dust mass-loss 
rates of the mass ejections that form the inner and outer dust shells 
are 1.8 x 107” M, yr ‘and 2.1 x 10° M, yr”, yielding total dust masses of 
3.2 x 10° M, and 5.2 x 10° M,,, respectively. 

However, newly formed dust is required to explain the substantial 
elevation of the MIR flux excess at day >+310. The fit results of the MIR 
excess of SN 2018evt are also shown as dotted red curves in Fig. 4. Our 
fit result is achieved by assuming the newly formed dust is composed of 


graphite grains of radius a = 0.3 um. In Fig. 6, we also present the mass 
of the newly formed dust as a function of time for a = 0.05 pm graphite 
and silicate grains. Dust masses estimated for other well-sampled 
supernovae that exhibit ejecta-CSM interactions are also presented for 
comparison, including the Type IIP supernovae 2004et (diamonds)"* 
and 1987A (blue upward triangles)"” and lIn supernovae 2005ip 
(crosses)”’”’, 2006jd (squares) and 2010jl (green downward trian- 
gles)”, As shownin Fig. 6, the amount of dust formed by SN 2018evt 
is equivalent to that formed in CC supernovae. 

We also remark that at day <+310, the MIR flux excess measured 
in band 2 (CH2 or W2) is higher than or comparable to that in band 
1 (CH1 or W1), indicating a higher dust emission efficiency towards 
longer wavelengths. This is compatible with the large (a =1.0 um) 
graphite dust particles in the primordial CSM shells suggested by our 
single-shell and double-shell models to the time evolution of the MIR 
excess at day <+310. 

Thea =0.3 um graphite dust model provides satisfactory fits to the 
MIR photometry in both bands 1 and 2 at day >+310 (Figs. 2 and 4). The 
indicated best-fitting radius of the newly formed dust grains also falls 
within the 0.01-1 um range of the typical size of the graphite dust grains 
(for example, refs. 107,108). However, the species of the newly formed 
dust grains may still not be inferred based on our observations as no spec- 
tralsignatures of CO overtone bands at 2.3-2.5 ppm were seen from our NIR 
spectra shownin Extended Data Fig. 3 (see also SN 2017eaw (ref. 125). Addi- 
tionally, we are not aware of any observation of SN 2018evt conducted at 
9 um, which may discriminate the silicate and graphite dust models", 
Therefore, we also present the results computed for a= 0.05 pm silicate 
and graphite dust grains in Fig. 6 and Extended Data Table 1. 


The progenitor’s mass loss 

Before the SN explosion, the progenitor mass-loss rate M can be associ- 
ated with the bolometric luminosity via a factor e, which denotes the 
kinetic-to-radiation energy conversion efficiency. Assuming a steady 
stellar wind CSM (s = 2in Pesm ~ F (ref. 29)), the bolometric luminosity 
Lpa can be written as 


v3 (3) 


where E,;, represents the kinetic energy of the thin shocked 
shell. The efficiency factor € is often assumed between 0.1 and 0.5 
(refs. 29,108,126,127). We adopted € = 0.3 and V,=2,000 kms“, the 
latter being consistent with the typical FWHM velocity of the inter- 
mediate Ha component measured over our spectroscopic campaign 
on SN 2018evt. The wind velocity blown from the progenitor was 
taken from the P Cygni feature reported in ref. 20, V, = 63 km s”. A 
similar velocity was observed only in the unshocked CSM of PTF11kx 
(ref. 128) (V, = 65 kms”), an 1999aa-like SN, which exhibits multiple 
CSM components but displays no signature of the ejecta—CSM interac- 
tions based on the early time observations’*"”’. 

A sudden decrease in the optical light curves of SN 2018evt can be 
seen at day ~+530, indicating the formation of new dust grains in the 
CDS (for example, Fig. 1 and Extended Data Fig. 4). Following the pre- 
scription in ref. 20, we approximate the optical bolometric luminosity 
(Lop) Of SN 2018evt by integrating its SED at day +530 over the optical 
wavelength range 3,870-9,000 A. The day +530 SED was obtained by 
warping the day +264 flux spectrum to match the BVgri-band pho- 
tometry at day +530. Therefore, the estimated Lop = 5.2 x 10" ergs “at 
day +530 yields a mass loss rate of 


: L v, 0.3\ /2,000 km s7! \? 
M x 0.04M, yr! bol )( w )( )( r ) 
oy (= x 10% ergs-!/\63kms-! € V; 


(4) 


The mass of shocked CSM around SN 2018evt can be estimated 
by multiplying the mass-loss rate to the duration of the shock 
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propagation (tjuration) aS approximated by the phase of the measurement 
tguration = 530 days, which can be expressed as 


Vee L 
M, = AMX tduration © 2.0M, (a ) 
shockedCSM V, duration O (52x104 ergs- 


( 22) (22% kms! j 52 ( fduration ) 
€ V; 530 days 


(5) 


At such late phases of SN 2018evt, the dominant radiation 
source inthe IR can be well-attributed to the thermal emission of newly 
formed dust (Fig. 2). Thanks to the MIR observations at day +517, a phase 
comparable to +530, we estimate the optical-to-MIR pseudobolometric 
luminosity Lope Of SN 2018evt by integrating the SED over a 
wavelength range of 3,870 A-5 um (Fig. 2). The computed 
Lopen = 1.2 x 10“ erg s™ at day +517 indicates a mass-loss rate 
M x 0.09 M, yr”. Therefore, the corresponding mass of the shocked 
CSM can be estimated to be M,nockea csm ~ 4-5 Mo- 

At day +517, adopting a shock velocity V, = 6,000 kms‘ estimated 
by the FWHM width of the broad Ha component (ref. 20 and Fig. 3) 
and an optical bolometric luminosity Lop: = 5.2 x 10” erg s”, following 
equations (4) and (5), the corresponding mass-loss rate and the mass 
of the shocked CSM yield M x 0.001 M, yr‘ and Mgpockeacsm = 0-2 Mo, 
respectively. If we include the emission in the MIR by adopting 
Loptemir = 1.2 x 10” erg s”, the corresponding M and Mgnockeacsm Can be 
estimated to be 0.003 M, yr ‘and 0.5 M,, respectively. 

In the literature, the luminosity of Ha line Lua serves as a good 
indicator of the bolometric luminosity L,,, as Lua has found to be pro- 
portional to Lpa (refs. 102-104,130,131). Lua can be expressed as 


Lita = 5 Ha Ve VR (6) 


where ep, denotes the efficiency of the conversion of the dissipated 
kinetic energy into H, luminosity in the shock wave. For SN 2018evt at 
day -+530, we measured Lya = 3.4 x 10*° erg s” (Extended Data Table 2). 
Thus we can get 


EH, _ Ly 
€ a 


Lua Lopt+MIR 
L z 0.03 x 
OPHEMIR (sa x 104° erg s-1 ) ( 1.2 x 1042ergs-! ) 


(7) 


We estimate Ea = 0.01 for SN 2018evt, which is comparable with 
the canonical value of e4, = 0.05 assumed in the literature’®”””®. 


Additionally, our computed a x 0.07 is also in general agreement 


Opt 


with that reported in Fig. 7 of ref. 20. 


The dust contributions of host galaxies by la-CSM SN events 
The far-infrared observations of the elliptical/lenticular (E/SO) gal- 
axies by the Herschel Space Observatory suggest that the typical 
dust mass of such galaxies spans 10*-10’ M, (refs. 132,133), while the 
average dust mass found in all types of galaxies in the local universe 
is 10°?" M,. Many dwarf elliptical galaxies exhibit dust masses less 
than 10° M, (ref. 134). 

Based on our optical-to-MIR observing campaign on SN 2018evt 
extended to day ~+1,000, we suggest that a total amount of -0.01 M, 
newly formed dust is formed in the postshocked region of the CDS. 
As the SN ejecta cools, more cold dust can be expected to form, as 
illustrated in Fig. 6. We remark that the 0.01 M, new dust formed in SN 
2018evt is estimated only for warm dust based on the 3.6 and 4.5 um 
observations by NEOWISE and Spitzer. If we assume it estimates the 
typical mass of the warm dust formed in la supernovae, the total mass 
of the newly formed dust could be higher by a factor of -10 if the bulk 
of the dust cools below ~30 K. Studies based on the observations of 
the Infrared Astronomical Satellite (IRAS) suggest ~90% of the dust in gal- 
axies was missed by IRAS, as IRAS and Spitzer are sensitive to warm dust". 


Depending on the detailed physical conditions, the timescale of 
the grain destruction could be as long asa few Gyr, based on the revised 
self-consistent models of dust destruction efficiency of supernovae’”° 
and other cases”. Moreover, the rate of la SN per unit mass decreases as 
the stellar mass of the galaxy increases (for example, Fig. 5inref.138).In 
particular, la SN rates of 6 x 10” and 6 x 10 * M, ‘year ‘are estimated 
for galaxy stellar masses of 10° and 10" Mg, respectively. Assuming 
0.1M, of cold dust (<30 K) is produced per Ia SN, la supernovae can 
produce onthe order of 10°-10° M, of dust for typical elliptical galaxies. 
Given the uncertainties in dust mass in ellipticals (M; < 10°-10’ M,), la 
supernovae can be responsible for 10% to 100% of all the dust in ellipti- 
cal galaxies. Considering that the la-CSM SN rate is about 0.02% to 0.2% 
of allla supernovae (ref. 18), the dust from la-CSM supernovae may be 
proportionally lower than the above estimate for supernovae la and 
cannot be the dominant source of dust in elliptical galaxies. We note 
that the effect of galaxy merging is also a dust source in E/SO galaxies 
as E/SO galaxies can capture younger galaxies together with their dust. 
The captured dust is usually distributed in a thin disk, but dust is also 
present ina diffuse environment (for example, ref. 139). The la-CSM SN 
contribution may also explain the diffuse dust. 


Data availability 

The data that support the findings of this study are openly available in 
the Science Data Bank at https://doi.org/10.57760/sciencedb.07968 
(ref. 140) or resolve.pid21.cn/31253.11.sciencedb.07968. The global 
network photometry at Las Cumbres Observatory is also available in 
the figshare repository https://doi.org/10.6084/m9.figshare.21543558. 
All spectra will also be made publicly available via Weizmann Interac- 
tive Supernova Data Repository. ATLAS*®” forced photometry service 
is available at fallingstar-data.com/forcedphot/queue/. The All-Sky 
Automated Survey for Supernovae* sky patrol interface is available 
at asas-sn.osu.edu/. The Spitzer Heritage Archive is available at irsa. 
ipac.caltech.edu/applications/Spitzer/SHA/. The NEOWISE co-added 
images and ALLWISE source catalogue are available at irsa.ipac.caltech. 
edu/applications/ICORE/ and irsa.ipac.caltech.edu/cgi-bin/Gator/ 
nph-scan?submit=Select&projshort=WISE. The Pan-STARRS data- 
base” is available at catalogs.mast.stsci.edu/panstarrs. Source data 
are provided with this paper. 


Code availability 

The SFFT” package used for image subtraction is publicly available 
at github.com/thomasvrussell/sfft. The IRAF software used for 
spectra reduction is available at iraf.net/. The BANZAI automatic 
pipeline used for image reductions at Las Cumbres Observatory is 
available at Ico.global/documentation/data/BANZAIpipeline/. The 
SExtractor® and PSFEx” software used for photometry are available 
at www.astromatic.net/software/sextractor/ and www.astromatic. 
net/software/psfex/. 
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Extended Data Fig. 1| The early-time comparisons of SNe 2018evt, and 1991T. rise time t, = 18.76 + 0.24 days, and a power-law index n =1.57 + 0.07 (grey curve). 


The spectrum of SN 2018evt (black curve) closely resembles the spectrum of SN The estimated t, is consistent with the V -b and rise time t,(V) = 20.00 + 0.68 days 
1991T at day -9 (red curve) in panel (a)”*”’, which exhibits the strong Fe IIIA 4404, of SN1991T/1999aa-like events”. The interpolations of BVR/light curves of SN 
A5129 absorptions and Ha emission, visible Si III 4564 absorption, Ni II blends 1991T are shown in dashed curves. The inset in panel (b) compares the B -V color 
around 4750A, 5350A, and weak S II W and Sill A 6355 as marked (for example,”*). curves between SNe 2018evt and 1991T, indicating a color difference < ~ 0.1 mag 
Panel (b) compares the early-time photometry of SN 2018evt (stars), and SN at similar phases. The corresponding Milkyway extinction is 0.05 mag™ for SN 
1991T (solid circles)”. A power law f ~ (t + t,)" is applied to fit the early V- and 2018evt. The error bars shown represent 1 - ø uncertainties of magnitudes, and 
c- band photometry, where r= (t-tg™)/(s(1+z))”?, tg "= 58352, s = 1.0 for colors. 


stretch value and z= 0.02523 for the redshift of SN 2018evt. The fitting yields a 
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Extended Data Fig. 2 | Optical spectra of SN 2018evt. Panel (a) shows optical the Ha profile of SN 2018evt from panel (a). For each epoch, the red dashed 
spectra of SN 2018evt spanning from days +136 to +579 relative to B-band line mirrors the spectral profile of the blue side across the peak flux of the 
maximum. Phases and facilities are marked on the right. Spectra of SN 2018evt intermediate Ha (for example, see Extended Data Fig. 6 for two Gaussian fits 
obtained at days +136, +272, and +351 are also compared to that of other Type to Ha). Its deviation from the red emission wing illustrates the time-variant 


la-CSM SNe (PTF11kx”*"” and SN 2002ic") at similar phases. Panel (b) portrays asymmetry of the Ha profile. 
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Extended Data Fig. 3 | NIR spectra of SN 2018evt. Phases and facilities are marked on the right spanning from ~ +125 to +516 days relative to the B-band maximum. 
Several most prominent lines are labeled. The near-IR spectrum of SN 2012ca obtained at days +198 - +208 is shown for comparison”. 


+516D Spex 
=23 
50 


Nature Astronomy 


Article https://doi.org/10.1038/s41550-024-02197-9 


eB-4 eg eit4 ae El 02ic ae hy aaa 11kx Ga Pee 13dn 
V E 2 ert 2 ire C ber 05g} au Aer 12ca 
a 
v -19 
Ke) 
14 2 
= 
D -18 
© 
= 
aa) 
16 < 7 
(ep) 
+ 
v = 
Ke) O 
E i Fs 
ea. <q ~°f © 18evt (CH1) 
@ 18evt (W1) 
= 
z C 
v w -19 
(49) O 
Q 20 = 
a = 
< c 
© 
-18 
z 
22 an) 
< 
(e9) 
+ 
5-17 
O 
e Las Cumbres “v 
24} v XingLong <x © 18evt (CH2) 
m Lijiang n @ 18evt (W2) 
100 250 400 550 700 100 300 1000 
Days since maximum light Days since maximum light 
Extended Data Fig. 4 | The optical and MIR light curves of SN 2018evt. Panel the light curves before and after day +310. The MIR light curves of several other 


(a) The BV gri band light curves of SN 2018evt. Panels (b) and (c) present the ~ 3.5 SNe la-CSM at similar phases are shown for comparison, including SNe 2002ic, 
pm (Spitzer CH1 and NEOWISE W1) and - 4.6 um (Spitzer CH2 and NEOWISE W2) 2005gj”, PTF11kx®’, 2012ca, 2013dn*"”, and 2020eyj”. The error bars shown 
photometry of SN 2018evt, respectively. Black-solid lines show polynomial fits to represent 1-ouncertainties of magnitudes. 
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NEOWISE. The first and the second rows display the Spitzer CH1 (3.6 zm) and images obtained from 2019-01-18 (day +149) to 2021-06-28 (day +1041) are shown 
CH2 (4.5 um)-band images obtained from 2019-05-19 (day +271) to 2019-11-09 in the remaining subpanels as labeled. Panels (a) and (b) display the Spitzer and 
(day +445), respectively. The third and the fourth rows present the NEOWISE W NEOWISE images, respectively. In each subpanel, the magenta cross indicates the 
1(3.4 um) and W2 (4.6 um) observations of the SN 2018evt field, respectively. location of SN 2018evt. North is up, east is to the left. 


The left column shows the reference images constructed by coadding the pre-SN 
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Extended Data Fig. 6 | The Ha profile of SN 2018evt observed with WiFeS at 
day +307 fitted with two Gaussian functions. The FWHM widths of the broad 
(cyan-dashed line) and intermediate (blue-dotted line) components yield 
5877 + 32 km s” and 1643 + 12 km s7, respectively. The red-solid curve gives 
the combination of these two components. The bottom gray line represents 
the Ha profile after subtracting the broad and intermediate components. An 
arbitrary offset has been applied to the residual spectrum for the purpose 


of presentation. The inset provides a zoom-in view of the P-Cygni profile as 
displayed in the residual spectrum. A double-Gaussian component fit to the 
residual spectrum near the Ha core is illustrated by the cyan curve. The location 
of the peak of the emission component suggests a redshift z =0.02561 + 0.00019, 
the location of the minimum of the absorption component measures a wind 
velocity V„=91 +58 kms”. 
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Extended Data Fig. 7 | Galactic extinction-corrected B-V, g-r, and g-i color are identical to the values at day ~ +310, as indicated by the horizontal black line 
curves of SN 2018evt. All colors were binned for 20 days to increase the signal- segments. The results of the 0.05 um graphite are shown by gray circles and 
to-noise ratio. The colors predicted by the scattering of the newly-formed linearly fitted by gray-dashed lines. The predicted colors of the 0.05 um silicate 
dust in the post-shock CDS are also presented by gray symbols as labeled. The dust are presented by gray squares and linearly fitted by gray-dotted lines. 
calculation was carried out by assuming that the intrinsic colors of SN 2018evt The error bars shown represent 1- g uncertainties of colors. 
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Extended Data Table 1| The MIR AB magnitudes of SN 2018evt, and the parameters describe the newly-formed dust 


0.3 um Graphite 0.05 um Graphite 0.05 um Silicate 
UT Date MJD Phase Map {CHi=3.6um) AB |CH2=4.5um) Temperature Dust mass Temperature Dust mass Temperature Dust mass 
(yy-mm-dd) (day) (mag) (mag) (K) (Mo) (K) (Ma) (K) (Mo) 


2019-05-19 58622.91 +270.91 18.20 + 0.08 17.92 +0.08 
2019-06-03 58637.50 +285.50 18.25 + 0.08 18.04 +0.08 
2019-10-29 58785.62  +433.62 17.50 +0.06 17.43 +0.06 997 +234 1.4+0.7x 10-4 1034+259 4.5+2.2x10-4 14664378 7.4+3.0x 1074 
2019-11-09 58796.79 +444.79 17.38 +0.06 17.34 +0.06 980 +77 2.0+0.4 x 1074 1026+89 6.1+1.4x1074 1511+193 9.142.2x 1074 


UT Date MJD Phase MAB.|W1=3.4um]  MAB,[W2=4.6um)] 
(yy-mm-dd) (day) (mag) (mag) 
2019-01-18 58501.21 +149.21 17.32 +0.09 17.33 +0.15 


2019-06-28 58662.37  +310.37 18.67 +0.20 18.12+0.21 
2020-01-20 58868.60  +516.60 17.07 + 0.05 16.89 + 0.07 706 + 39 1.740.4 x 10-7 727443 5.541.2x 10-3 918+68 9.4+2.0x 10-3 


2020-06-26 59026.36  +674.36 16.81 +0.05 16.63 + 0.07 689432 2.8+40.5x10-3 707434 9.0+1.7x10-3 892458 1.5+0.3x 10-7 
2021-01-18 59232.70 4880.70 17.53 +0.06 16.93 +0.10 5544+43  6.6+1.3x107?  566+43  2.10.8x107? 675+43 = 3.64 0.8 x 107? 
2021-06-28 59393.48 +1,041.48 18.43 +0.11 17.43 +0.08 467433  1.2+0.2x107? 4764+33 3.7+1.6x107? 5494+33 6.741.6x 107? 


The CH1 and CH2 magnitudes were measured from the images taken by the Infrared Array Camera (without host subtraction) and the host-subtracted W 1 and W2 magnitudes were measured 
from the observations by the NEOWISE reactivation mission. The presented parameters are deduced for a = 0.3 um, and a = 0.05 um graphite and silicate dust grains. 
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Extended Data Table 2 | Log of the optical and NIR spectroscopic observations of SN 2018evt 


UT Date! MJD Phase? Resolution Range Instrument/Telescope Exposure time Airmass logliyg EWea 11 TIR triplet 
(yy-mm-dd) (day) (Å) (Å) (s) (log(erg s™!)) (A) 
2018-08-12* 58343.00  -—9.00 15.8 3650-9200 EFOSC2/NTT 300 1.42 a 2 
2018-12-24* 58476.64 124.64 ~ 15.0 3500-10000 FLOYDS/2.0-m FTN 1800 1.49 41.53+0.01 —846.90 + 5.01 
2019-01-01* 58484.70 132.70 ~ 15.0 3500-10000 FLOYDS/2.0-m FTS 1600 2.09 41.50 +0.01 —775.80 +4.75 
2019-01-04 58487.94 135.94 25.0 3400-9100 YFOSC/LJT 1350 1.38 41.48 +0.01 ass 
2019-O1-11* 58494.58 142.58 ~ 15.0 3500-10000 FLOYDS/2.0-m FTN 1600 1.56 41.52 +0.01 —774.40 + 4.64 
2019-01-21* 58504.64 152.64 ~ 15.0 3500-10000 FLOYDS/2.0-m FTN 1600 1.18 41.52 +0.01 —857.92+5.11 
2019-01-31 58514.89 162.89 25.0 3400-9100 YFOSC/LJT 1350 1.29 41.52 +0.01 WA 
2019-03-04* 58546.43 194.43 ~ 15.0 3500-10000 FLOYDS/2.0-m FTN 1800 1.62 41.55 +0.01 —1093.59 6.55 
2019-02-08 58550.75 198.75 15.0 4000-9000 BFOSC/XLT 3300 1.92 41.45 +0.01 saa 
2019-03-17* 58559.49 207.48 ~ 15.0 3500-9250 FLOYDS/2.0-m FTN 1800 1.18 41.54 +0.01 — 1082.88 + 6.53 
2019-02-20 58562.75 210.75 15.0 4000-9000 BFOSC/XLT 3600 2.46 41.46 +0.01 

2019-03-28 58570.73 218.73 25.0 3400-9100 YFOSC/LJT 1500 1.31 41.50+0.01 sae 
2019-03-30* 58572.49 220.49 ~ 15.0 3500-10000 FLOYDS/2.0-m FTN 1800 1.17 41.51 +0.01 —1151.35 +6.96 
2019-04-13 58586.75 234.75 15.0 4000-9000 BFOSC/XLT 2700 1.63 41.47+0.01 

2019-04-23* 58596.66 244.66 ~ 15.0 3500-10000 FLOYDS/2.0-m FTS 2700 1.31 41.49+0.01 as 
2019-04-26 58599.30 247.30 15.8 3650-9200 EFOSC/NTT 1500 1.51 41.47+0.02 —1302.33+ 14.06 
2019-05-11* 58614.39 262.39 ~ 15.0 3500-10000 FLOYDS/2.0-m FTN 2699 1.19 41.44+0.01 —1257.79 + 7.65 
2019-05-13 58616.22 264.22 15.8 3650-9200 EFOSC/NTT 1499 1.28 41.43+0.02 —1545.90+ 16.17 
2019-05-20 58623.69 271.69 25.0 3400-9100 YFOSC/LJT 2000 1.34 41.40 +0.01 ast 
2019-06-09 58643.14 291.14 15.8 3650-9200 EFOSC/NTT 1800 1.24 41.38+0.02 —1387.15 + 14.73 
2019-06-09* 58643.37 291.37 ~ 15.0 3500-10000 FLOYDS/2.0-m FTN 2700 1.40 41.39+0.01 —1162.68 + 7.25 
2019-06-24 58658.50 306.50 ais 3000-9500 WiFes/ANU 1200 es 41.47+0.01 as 
2019-07-15* 58679.25 327.25 ~ 15.0 3500-10000 FLOYDS/2.0-m FTN 3600 1.36 41.30 +0.01 —1221.3447.56 
2019-08-08 58703.04 351.04 15.8 3650-9200 EFOSC/NTT 2699 1.88 41.26+0.02 —1296.33 + 14.30 
2019-08-22* 58717.37 365.37 ~ 15.0 3800-10000 © FLOYDS/2.0-m FTN 3600 1.76 41.22+0.01 —1145.15+7.78 
2019-12-24 58841.62 489.62 ~ 15.0 3500-10000 | FLOYDS/2.0-m FTN 3600 1.54 40.76 +0.01 — 1036.41 + 6.83 
2020-01-19 58867.57 515.57 ~ 15.0 3500-10000 FLOYDS/2.0-m FTN 3600 1.42 40.61 +0.01 —877.62 +6.46 
2020-02-03 58882.54 530.54 ~ 15.0 3500-10000 FLOYDS/2.0-m FTN 3600 1.34 40.53 +0.01 —1071.95 +7.56 
2020-02-19 58898.29 546.29 15.8 3650-9200 EFOSC/NTT 2700 1.12 40.37 +0.02 

2020-03-23 58931.31 579.32 15.8 3650-9200 EFOSC2/NTT 2700 1:12 40.16 + 0.02 j 
2018-12-24 58476.65 +124.65 ~ 18.3 7000-25000 SpeX/IRTF 150 x 10 1.45 wi —586.80 + 11.99 
2019-01-01 58484.56 +132.56 ~ 16.0 8000-25000 | GNIRS/Gemini North 90 x 20 1.69 ar —752.82 + 5.06 
2019-01-08 58491.63 +139.63 ~ 16.0 8000-25000 GNIRS/Gemini North 90 x 20 1.28 os —803.16+5.53 
2019-05-11 58614.40  +262.40 ~ 18.3 7000-25000 SpeX/IRTF 150 x 10 1.17 re —1051.25 + 22.74 
2019-05-16 58619.22 +267.22 ~ 24.0 7800-25000 FIRE/Magellan Baade 126.8x 12 1.26 zs — 1136.63 +8.47 
2019-06-17 58651.30 +299.30 ~ 18.3 7000-25000 SpeX/IRTF 150 x 10 1.17 sa —1107.76 + 23.33 
2019-07-12 58676.03 +324.03 ~ 24.0 7800-25000 FIRE/Magellan Baade 126.8x8 1.13 oe —1206.97 + 8.58 
2020-01-19 58867.64 +515.64 ~ 18.3 7000-25000 SpeX/IRTF 150 x 10 1.18 si —761.89 + 17.06 


The measured Ha luminosity and the equivalent width of the Ca II NIR triplet are also listed. a* marks the spectra that are already published in*°*. b Days since B-band maximum on MJD 58352 
/ 2018 August 22. c Uncertainty is derived and assumed that all spectra have 10% flux uncertainty. Note that the distance is not included. Only the Milky Way extinction is corrected with E(B-V) 
MW =0.05 mag. R= 311. 
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